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This paper proposes a decentralized constrained optimal control of the 
multimachine power system stability. Today’s power networks operate closer 
to their maximum limits. Alternative Current power grids are more vulnerable 
and subject to instability than ever before. A three machine power system and 
four machines, power system connected with a lossy transmission line, is 


investigated. The linearized dynamical equations of the multimachine power 
system are obtained near to an equilibrium point, and it can be stabilized by 
using a decentralized constrained controller based on optimal control. The 
feedback controller, which comprises independent control stations, receives 
the measurement data and influences the control input of the machine attached 
to it the subsystems. State feedback controller guarantees that the closed-loop 
system is stable and can guarantee the performance index. The designed 
controlled systems are based on the Algebraic Riccati Equations, and all its 
poles are in the closed left half-plane. Decentralized constrained optimal 
control of the multimachine power system is achieved through simulation 
results. The following results are proposed by improving power system 
stability. 
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1. INTRODUCTION 

Power systems stability problem has been enumerated and continues to receive a great deal of attention 
over the years. Stability of a power system is the ability of the power system to return to its normal conditions 
after being subjected to a disturbance, such as a system fault such as a transmission line outage for generator 
outage [1-2]. Supporting a power system with several machines and supplying continuously that is a motivating 
operating system in different scenarios could not lead to failure of the procedures and blackouts. A 
decentralized controller is a feedback controller, which consists of independent control stations, each of which 
receives the measurement data and influence the control input only from the attached subsystems [3]. The 
information flow from the plant through the controller in the feedback loop that is divided into a separated 
parallel flow through the control stations. Decentralization is an effective way to deal with complexity and 
cross changes in the interactions of large-scale dynamic systems. Although the optimal control concept has 
many advantages, it involves a complicated route in solving the problems of nonlinear systems. 

Different methods have been used to improve decentralization of the multimachine in recent decades. 
Xi [4] designed a nonlinear decentralized controller used in Hamiltonian functions. The methods of the 
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saturated steam valve and excitation controller in the power system. Ortega [5] Proposed transient stabilization 
of multimachine power systems with nontrivial transfer conductance based on Hamiltonian port and estimated 
the domain of attraction of a bona fide Lyapunov function. Li [6] proposed a nonlinear robust decentralized by 
using linear matrix inequalities approach to multimachine cases, based on suboptimal guaranteed stability, cost 
performance to be operated by different operating point fault location. Abouelsoud [7] proposed a stabilizing 
controller for a single machine power system with a nonzero conductance of transmission line, based on the 
observer to estimate the power angle of the system. The stability of the reduced-order observer and closed-loop 
of the system incorporating using the Lyapunov direct method is for a single machine. 

Huang [8] designed a new model of excitation decentralized constrained controller of the 
multimachine power system based on the optimal feedback gain. Lee and Liaw [9] proposed the solution of 
inverse linear optimal control for constant disturbance using Algebraic Riccati Equations based on the linear- 
quadratic regulator for an optimal feedback law and the relationships between the open-loop poles and the 
closed-loop poles that guarantee a positive regulator. Abouelsoud [10] proposed a global stabilization of linear 
time-invariant systems subject to state and control constraints based on the direct feedback controller to 
globally asymptotically stabilize a continuous-time linear time-invariant subjected to symmetrical polyhedral 
state and control constraints sets. 

Elloumi [11] compared a decentralized control law to implement a multimachine power system. 
And it presented two aspects which are a linear decentralized control, and the gains depend on nonlinearity of 
the system and successive approximation approach to determine two-point boundary value. Hongshan and 
Xiaoming [12] proposed a fast excitation predictive control method for a multimachine power system based 
on a dynamic model and some inequality constraints on states, inputs, and outputs are considered in rolling 
optimization by using Gramian balanced reduction technique. Hosni and Selwa [13] investigated the problem 
of decentralized robust stabilizing control, based on direct Lyapunov methods and Kronecker reduction 
production notation. Tao [14] proposed an optimal control problem for large-scale systems with unknown 
parameters and dynamics. It uses a robust adaptive dynamic programming method and decentralized optimal 
control with unmatched uncertainties. The convergences of the robust adaptive dynamic programming 
algorithm and asymptotic stability of the closed-loop of the large-scale system are studied rigorously. 

This paper proposed an improvement in the multimachine power system using decentralized 
constrained optimal control. Transient stability assesses using parameter space has been limited by the number 
of machines and stability are bounded because of size and dimension. Power system stabilizer and genetic 
algorithm have applied to a single machine infinite bus, which complicated for the multimachine system. The 
difficulties arise that use a random variation and selection of minimizing estimate error and time consumed for 
the artificial immune system and evolutionary programming [23-25]. 


2. MODEL AND PROBLEM FORMULATION 

The paper considers the multimachine power systems consisting of generators interconnected through 
a transmission line which is assumed as a lossy. The dynamics of the i“ generators is presented 
by [5, 15-16]. 

The classical three-dimensional flux decay model presents the dynamic of the i" generator with 
excitation. 








5-0, (1) 
* — Dm , Do 

=a tg (Pa Pa) 

E=—(-E,tE,tvieN= th siueseiice N} 


Where the states variables of this subsystem are the rotor angle 6i€ I: =[0,272), rad, the speed deviation wi € 

R in arad/sec and the generator quadrature internal voltage Ei € R: in p.u. Dp is the damping constant in pu. 

Mi= 2H, where H is the inertia constants in seconds, Pmiis the constants mechanical input power and P,i is the 

electromagnetic power, both in pu (per units), t; is the direct axis transient short circuit time constant. The 

voltages Eqi, Es and vi, expressed in pu are the electromotive force in the quadrature axis, 

the control constant voltage component applied to the field winding and the control voltage input, respectively. 
The active Pj and reactive Qei power, as well as the voltage Eqi, are defined as: 
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(2) 


Where I, and Ig; are the quadratures and direct axis currents, Xqi is the direct axis reactance and x’ qi the direct 
axis transient reactance, where it says that Xai>x’ ai. Ig is the field winding current and x,qi 1s the mutual reactance 
between the excitation and stator coils. All the latter quantities are expressed in p.u. System. 

The currents give the interconnection between machines. 


16" LE (Gs06,*BasS"5;) (3) 
1s LE (Gui885,~ Buse) 


Where 63; = 6;-6;, Gmij is the conductance and Bnijis the susceptance both in pu resulting from the computation 
of the network admittance matrix. Developing the sums when j=i and recalling the identities. 
The network admittance through conductance and susceptance reads: 


Y, = GiB 
Gas (4) 


mij 


a ii = tn 


GniiSi®. 5 Bry Oy = —Y 4©5(5 + ats) 


Yields the current in quadratic and direct axis: 
N’ . 
1,=GuE* DEY ,Si(5,+a1,) (5) 
iz 
N 
Le = -BwEi- EY (5 ;+a;) 
j=l 


Finally, combining (1) and (4) results in the well-known compact form is the dynamic equation of the 
third-order model (5 1» @,, E. ) which becomes: 


01-0; 
. ‘ N ; 
QQ; Da; P, GF: dE, > EY, sin(5,+a,) 


j=l, ji 


(6) 
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With the positive constants: 





D.. 
D, = Pp aa dP, .5 Gi; = dG i> = ae 
M., M 


i i 


1 ; 1 ' 
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Similarly, it observes that, ai, bi>O aj=ajj, and that if M=Mj, Yi=Yji 
It assumes that the model with u;=0 has stabled the equilibrium points [6..,0, E, | and it finds the control law uj 


such that it operates in closed-loop. 
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3. STABILIZATION WITH STATE CONSTRAINTS 
Considering a continuous-time Linear Time-Invariant system as follows: 


x(t) = Ax(r)+Bu(t) (7) 


x € R",u € R™, (A,B) it is a controllable pair and symmetric constraint state and control of sets. 





S, = {x € R":-a, < Gx < w,} (8) 











S, = {u € R":-a, < G,x < w,} (9) 


G, € R°1*", E,, € R™*" Are both full row rank, it considers the problem of designing a state feedback 
controller [10, 17, 21, 22]. 

By scaling, let us make w, = 1 andw, = 1 

S, = {x € R":-1<G4,x < T} 

$= eR 1 =F x= 1} 


=[&] -a(° ) 
c=[@)2= Ce 
u = K,, then the state and control constraints become 
<2 (s,)#S10r-1<6, +8, <1 
Assume that the invariant zeros of the system (A, B, G, and E) are in the open left-half plane. 

















u(t) =—Kx(t) (10) 


K Is gain, such that the closed-loop system becomes: 


x(1)=A.2x(1) (11) 


Where A, = A — BK is asymptotically stable, and both state control constraints (8), (9) are satisfied. This 
implies that the state Sx is A, invariant. 


3.1. Simulation Results 
The motion equation of the three-machine power system can be described by: 


6, =Q, 


@, =—D,o,+ P. —G,,E; —YE,E, sin(6,-6,+@) 


E, =—a,E, + b,E, cos(6,-6,+a@)+E,, +4, 

: (12) 
0, =Q, 

@, =—D,a, + P, —G,,E; — YEE, sin(6, — 6, —@) 
E, =—a,E, +b,E, cos(6, —6,+@)+E,.+u, 


3 = WD; 


@, =—D,o, + P,— G,,E; — YE,E, sin(6, — 6, — a) 


he 


, =-a,E,+b,E, cos(6,—-6,+a@)+E,,+u, 


During the steady-state, all these parameters become as mentioned below to determine the equilibrium points: 


o, =0.a, =0,a, =0.u,=0,u, =0.u,=0.a@=0. 


Oi Onis and 1 oie ; E., oer to be known and then; 
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(13) 
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2: . 
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Considering a power system network represented by four machines connected to a transmission line 
is lossy, as shown: 
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6, =Q, 


o, =—-D,o,+ P.—G,,E; — YEE, sin(6,—6, +@) 
E, =-4,E, +b,E, cos(6,—-6, +a) + E,, +, 
6, =, 


@, =—D,@, + P,—G,,E; — YE,E, sin(6,—6, -@) 





E, =—a,E, +b,E,cos(6,—6,+a@)+E,,+u, 


6; =O, 


o, =—D,a, + P,—G,,E; — YE,E, sin(6,—6,;—@) 





E, =—a,E, +b,E, cos(6,—-6,+@)+E,,+u, 


6, =, 


o, =—D,0,+ P,—-G,,E; —-YE,E, sin(6,—6,+@) (24) 





E, =—a,E,+b,E, cos(6,—-6,+@)+E,4+u, 


It determines the equilibrium points by assuming: w, =0; @, =0;@,=0 ; @,=0 ; u,=0; u,=0 ; u,=0 


;u,=0; The calculation of the parameters is: Di s0%y : Os Ow CE ties Fay EG 







0,010 + j0.085] 0.032 +j0.161 


0.017 +j0.092 


Figure 1. Three machines power systems [1] 
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Figure 2. Four machines power systems [18-20] 
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4. RESULTS AND DISCUSSION 

These results are compared and discussed in the numerical procedure for the multimachine power 
system. Decentralized constrained controller not only stabilizes the power system but also achieves an optimal 
solution for overall subsystems. The numerical simulation of the three-machine power system and the four- 
machine power system indicated the results. 


4.1. Three Machine System 

The simulations carried out on the three machines, a power system with a lossy transmission line and 
decentralized constrained controller is solved. A three-phase fault occurs at 1.5 seconds on the transmission 
line between bus 7 and 8 nearest to bus 8 in Figure 1. The fault is cleared by opening the circuit breakers both 
at the ends of the line. It makes over 2 seconds. During the steady-state, the conventional control system with: 


Uy = 0 Uy = 0:u, =() , which cannot maintain the stability of the fault is cleared in 2.2 seconds after fault 


occurrences. 


The optimal control decentralized constrained, and they give the feedback state gain matrix K time 


response of the rotor angle; the rotor speed and the excitation system of the machines O, : O. : O. ;, ; O,;O; 


: E 1% Ez : E, are illustrated in the Figures 3 in Figure 11 which show the responses of the multimachine system. 


Those figures show the advantage of the proposed control scheme. The simulation of the proposed control 
structure of the system under identical conditions leads to the following observations. It plots the results that 
show the rotor angle versus time, the rotor speed versus time, and the excitation system versus time are shown 
in Figure 3 to Figure 11 for each machine systems. 

This control technique involves several independent local controllers decentralized through subsystem 
coordination. It considers a local feedback controller communicating among subsystem with constraint 
bounded. This advantage of the decentralized controller is reliable and enlargeable for each subsystem seems 
to achieve its objective. The decision is taken locally, and negotiation can take place in different actuators. 
They usually have different goals that are applied for three machine system and four machine system. 
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Figure 11. Excitation system of machine #3 


The analysis assumes that when there are no losses, the conductance is null. The equilibrium points at 


the steady-state after the fault is: (5,,.a,, £,.,5).,0),E>+5ysQysEy)= (0.03965,0,1.0566,0.344, 0, 1.0502, 0.2299, 0,1.0170) 
which verifies the conditions of the three machine power system in the dynamical model. The rotor angle of 
all the machines shown in the Figures 3, 6 and 9, are increasing during the fault occurrence and it regains to 
the normal state after clearing the fault between 0.7 seconds to 0.9 seconds, and it converges to the normal 
state. The rotor speed of the system in the Figures 4, 7 and 11, is in the steady-state and the increase after 
introducing faults at 0.7 seconds to 1 second and it regains to the normal state after the fault occurred. Voltage 


increased till 1 second and restore to the maximum level after the fault is cleared at 1 second to the normal 
state. 


4.2. Four Machines System 

This subsection presents simulations of the four-machine system in Figure 2. In this case, 
the disturbance is a three-phase bolted fault are applied in 0.5 seconds in the bus 9. This fault is cleared by 
tripping line 1. The evolution of the rotor angle, the rotor speed, and the excitation system of the generator are 
observed and detailed. Finally, graphs are shaped and presented. 





Decentralized constrained optimal control of the multimachine ... (Djibrine. Abakar) 


1166 O ISSN: 2502-4752 


That modifies the topology of the network and consequently induces changes in the equilibrium 


points, as identified: (5,0, , E,,,5).. Ores Ere Sjes Osos Ese, Ojos Oyes Ese) = (0.61085,0,1.0397,0.8039,0,1.16,0.6105, 0,1.0397, 0.8039, 0,1.16) 
and it verifies the conditions of the four machines power systems as connected. 

The rotor angle of the machines is presented in Figures 12, 15, 18 and 21, are stable during the steady- 
state and increase during the fault and stabilized at one second. The rotor speed of the machines is shown in 
Figures 13, 16, 19 and 22, are increasing up to 0.7 seconds, and it sustained at one second during the fault. The 
excitation system of the machines is depicted in the Figures 14, 17, 20 and 23, are decreased during from 0 
seconds to -0.5 seconds and to stabilize during 10 milliseconds at overall and so far, to regain at one second to 
the normal state. 
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Figure 16. The rotor speed of machine #2 Figure 17. Excitation system of machine #2 
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Figure 18. Rotor angle of machine #3 
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Figure 20. Excitation system of machine #3 
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Figure 22. The rotor speed of machine #4 


4.3. Analysis of the Results 
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Figure 19. The rotor speed of machine #3 
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Figure 21. Rotor angle of machine #4 








Figure 23. Excitation system of machine #4 


This section analyzes the performance of the closed-loop system. Which are subjects, both to a three- 
phase fault occurs at the transmission line between buses. Simulation of the proposed control structure of the 
multimachine power system, under this condition, leads to the following observations. It presents the resulting 
rotor angles versus times; rotor speeds versus times and the excitation system versus times are shown in the 


figures with corresponding. 


The control technique requires feedback from each subsystem. Decentralized through subsystem 
coordinating a local feedback controller is considerate to communicate among subsystem with constraints 
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bounded. The closed-loop system is asymptotically stable with approved performance measures. 
The advantage of a decentralized controller is reliable and expandable for each subsystem seems to achieve its 
objective. The decision takes place locally, and negotiation can take place in different actuators. 
The system exchanges and subproblems can usually be cooperative or competitive. 


5. CONCLUSION 

This paper presents decentralized constrained optimal control to enhance the stability of the 
multimachine power system. The main contribution and simulation results based on the Algebraic Riccati 
Equations using the decentralized constrained controller to improve the stability of the multimachine power 
system has been achieved. Different approaches explore it and compute the following conclusion: 

First, the stability of the multimachine power systems for state constraints bounded and all its poles are in the 
closed left-half plane. Second, the decentralized optimal controller has been proposed, and the feedback from 
the subsystems is to control themselves. 

Improving the stability of the multimachine power systems assumed a lossy transmission line to 
measure and estimate the values of excitation and angles. Designed optimal controller decentralized with state 
constraints bounded. Decentralization is an effective way to handle complexity and cross changes in the 
interaction of large-scale dynamic systems. It often controls the large-scale system by over one controller or 
decision-maker encompassing decentralized computations. 

It satisfies the stability of the overall closed-loop. In any control design, it pursues a controller not to 
stabilize the system outbuilding also to guarantee the adequate performance of the system. All the subsystems 
level to secure the entire subsystem asymptotically stable with a performance measure. It has shown that the 
controller is not aggressive. The paper has also demonstrated that the proposed results yield necessary and 
sufficient application to the optimal constraint control of three machines and four machine power systems and 
validate the findings. It has not only stabilized the multimachine power system but also completes the 
suboptimal control guaranteed, cost performance index for all bounds controlled. 
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